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Abstract 

We investigate the transmission and reflection of Bose-condensate excitations in the low energy 
limit across a potential barrier separating two condensates with different densities. The Bogoliubov 
excitation in the low energy limit has the incident angle where the perfect transmission occurs. 
This condition corresponds to the Brewster's law for the electromagnetic wave. The total internal 
reflection of the Bogoliubov excitation is found to occur at a large incident angle in the low energy 
limit. The anomalous tunneling named by Kagan et al. [Yu. Kagan et al, Phys. Rev. Lett., 90, 
130402 (2003)] can be understood in terms of the impedance matching. In the case of the normal 
incidence, comparison with the results in Tomonaga-Luttinger liquids is made. 

PACS numbers: 03.75.Kk, 03.75.Lm 
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I. INTRODUCTION 



The transmission and the reflection are basic concepts in physics, that have been stud- 
ied for strings, elastic bodies, fluids, electromagnetic waves, quantum particles, Tomonaga- 
Luttinger liquids and so on. The present paper is devoted to the study of the transmission 
and the reflection of collective excitations in the condensed Bose system. This proceeds 
from the earlier study of the so-called anomalous tunneling |l( , which is described as follows: 
Bogoliubov excitations transmit perfectly through the potential barrier separating two con- 
densates in the low energy limit. This is in markedly contrast to single particle tunneling in 
the usual quantum mechanics, where the perfect reflection occurs in the low energy limit. 

We review a rough sequence of studies of the anomalous tunneling. Kovrizhin initially 
discussed the transmission coefficient of Bogoliubov excitations through the 5-function po- 
tential barrier [2J|. Kagan et al. also studied the tunneling of Bogoliubov excitations through 

n n 

the rectangular potential Through the study in Ref. Kagan et al. pointed out that 
this phenomenon is anomalous, in comparison to the usual quantum tunneling of a single 
particle. 

Danshita et al. considered the transmission of Bogoliubov excitations across the 8- 
function barrier separating two condensates with different macroscopic phases |3|. They 
found that the anomalous tunneling disappears when the phase difference reaches the criti- 
cal value giving the maximum supercurrent of the condensate. 

These works were done using the one-dimensional potential. On the other hand, the 
problem of Bogoliubov excitations scattered by a spherical potential was examined 0, 3, E| • 
The energy dependence of the cross section of the low energy scattering by the spherical 
potential is consistent with the Rayleigh scattering of a sound wave in the classical wave 
mechanics, i.e., the cross section vanishes in the low energy limit 

Using the fact that the wave function of the excited state in the low energy limit cor- 
responds to the macroscopic wave function of the condensate [7], Kato et al. concluded 
that the anomalous tunneling occurs for a potential barrier V(x) being arbitrary symmetric 
faction („(-*) = VV) g They als0 Moated that the perfect tuning appears even 
at finite temperatures within the scheme of the Popov approximation [8J. 

In the Tomonaga-Luttinger liquids, on the other hand, the perfect transmission has been 
also discussed in the context of quantum wire, using the renormalization-group B], and 
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conning the mu.tip.e renews 10], EspeciaHy, using the renormahzation-group theotv, 
it was shown that the barrier is an irrelevant perturbation in an attractive fermion system |9[ . 
That is to say, the perfect transmission is known to occur in Tomonaga-Luttinger liquids in 
superconducting phase. 

As seen in the above, the property of the transmission of the Bogoliubov excitation in 
the low energy limit has not only a specific character of the condensed Bose system but also 
something in common with that of the classical wave mechanics and that of the Tomonaga- 
Luttinger liquids. From this point of view, we have three aims in this paper as follows: (i) 
exposing novel phenomena of the transmission and the reflection of the Bogoliubov phonon, 
(ii) giving a physical implication about the anomalous tunneling analogous to the classi- 
cal wave mechanics, (iii) comparing the result of the tunneling problem of the Bogoliubov 
phonon with that of the Tomonaga-Luttinger liquids. 

To achieve these aims, we investigate the tunneling problem of the Bogoliubov phonon 
through the arbitrary potential barrier separating two condensates with different densities. 
We consider that the incident Bogoliubov phonon does not only run toward the wall per- 
pendicularly but also runs toward the wall with an arbitrary incident angle. This extended 
problem is analogous to the problem of the reflection and the refraction of a wave at the 
interface between two different mediums. 

We may summarize the results of this paper as follows: (i) The Bogoliubov phonon has 
an incident angle where the perfect transmission occurs between condensates with different 
densities. This condition corresponds to the Brewster's law, i.e., the sum of the incident 
angle and the refracted angle is equal to n/2. There also exists the total internal reflection 
of the Bogoliubov excitation in the low energy limit, (ii) The perfect transmission of the 
Bogoliubov excitation in the low energy limit can be regarded as a result of the impedance 
matching between equivalent condensates separated by the potential wall. The impedance 
is inversely proportional to the sound speed of the Bogoliubov phonon. (iii) At normal inci- 
dence, our result in the low energy limit is consistent with the result obtained by the theory 
of Tomonaga-Luttinger liquids 10], when one uses the Luttinger liquid parameter for the 
weakly interacting Bose gases. We note that the negative density reflection in the interacting 
condensed Bose system cannot be necessarily identified with the Andreev reflection. 

The outline of this paper is as follows. In Section II, we shall give a formulation of the 
problem, and derive the transmission and reflection coefficients. To obtain these coefficients, 
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we use properties that the wave function of the excited state in the low energy limit corre- 
sponds to the condensate wave function, and the constancy of the energy flux. In Section 
III, we shall discuss results in accordance with three aims. We also propose future problems 
on the experimental and theoretical sides. In Section IV, we summarize our results. 



II. FORMULATION AND RESULTS 



In this section, we derive the transmission and reflection coefficients of the Bogoliubov 
excitation in the low energy limit. We use the mean-field theory, say, the Gross-Pitaevskii 
equation 
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12] and the Bogoliubov equation 13]. To discuss the property of the Bogoliubov 



excitation in the low energy limit, we use the fact that the wave function of the excited state 
in the low energy limit corresponds to the macroscopic wave function of the condensate [?J ■ 
To evaluate the transmission and reflection coefficients, we use the constancy of the energy 
flux. 

The stationary Gross -Pitaevskii equation written in dimensionless form is given by 

W*(r) = 0, (1) 

where 

7^-iv 2 + F(r)-l + |^(r)| 2 , (2) 

and the Bogoliubov equation also written in dimensionless form is given by 

h! -\Kf) 2 \ / u(f,e) \ / u(f;e) 




2 . || _ II _ I ' 

* (f) —h' I \ v(r,e) I \ v(r;s) I 
where h! = TC + ^(f)) 2 . We have introduced the following notations r = r/£, V = £V 



\l/(r) = g / fi^f (r) , V(f) = V(r)/fi, and e = s/fie, where g(> 0) is the coupling constant 
of the two-body short-range interaction, [i is the chemical potential, and £ is the healing 
length defined by £ = h/y/mjl with m being the mass. Henceforth, we omit the bar for 
simplicity. We take the potential V(r) to be a function of x and allow V(x) to be asymmetric 
V(x) 7^ V(—x). We consider V(x) being superposition of short-range potential near x = 
and a potential step with asymptotic form: 

V(x) <! VL ' f ° r X <K _1 (4) 
V n , for x » 1, 




FIG. 1: A schematic figure of the system. In this figure, we take xy plane as an incident plane. 

where both Vl and Vr are smaller than unity. The profile of V(x) is shown in Fig. [TJ 

There are steady-states of a Bose-Einstein condensate even in the presence of a potential 



step In this paper, we shall treat a case without a supercurrent where the phase is 

spatially constant. Henceforth, \I/(r) is considered to be real. In this situation, we have the 
following asymptotic behavior of the condensate wave function: 



VI - VL, for x < -1 

- { , (5) 

VI - V R , for x > 1. 

We note that our objective is to derive the transmission and reflection coefficients which can 
be adopted to the general case. However, we shall give an example here, in order to compare 
our general result with a concrete example. For instance, we assume that the potential V(x) 
has the following behavior: 



0, for x < 0, 



V(x) 



V b , for < x < 1/V2, (6) 
0.8, for 1/V2 < x. 



We use the strength of the potential barrier Vb as Vb = 3 and 5. In Fig. [2l we show 
condensate wave functions obtained by numerical calculations, in the presence of the above 
potential V(x). A solid line represents the numerical result for Vb = 3, and a dotted line 
represents the numerical result for Vb = 5. These numerical solutions satisfy the asymptotic 
behavior given in Eq. (jSJ). When the potential energy in the asymptotic regime is less than 
the chemical potential, i.e., V u < \i with v = L for x <C — 1 and v = R for x ^> 1, there are 
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FIG. 2: The numerical solutions of the condensate wave function. We use the potential energy 
as the function given in Eq. (|6|). A solid line represents the case for = 3, and a dotted line 
represents the case for Vb = 5. 

condensate wave functions spatially constant in the asymptotic regime, as shown in Fig. [2J 
Hence, Bogoliubov excitations even in the low energy limit exist in both sides of the potential 
barrier. 

We shall consider the situation where a Bogoliubov excitation with a wavenumber vector 
Icl runs against the potential wall at an angle 0l with respect to the normal vector of the 
wall (i.e., rc-direction). The Bogoliubov phonon is split into a transmitted wave with a 
refraction having a wavenumber vector at an angle 0r and a reflected wave with at 
an angle <ft' L as shown in Fig. [TJ 

Owing to the translational invariance in y, z directions, the Gross-Pitaevskii equation 
reduces to 




(7) 



with 




(8) 



The solution of the Bogoliubov equation has the form of 
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and the Bogoliubov equation is given by 

h" -^ 2 {x^ ( 
m 2 (x) -h" 



u[x\ e) 



v(x; e) 



(10) 



u(x; e) 
y v(x; e 

where h" = h + ^ 2 {x) + k\/2, and k± = y/K* + k\. 

Let us consider asymptotic forms of a Bogoliubov mode. At \x\ ^> 1 where ty(x) 
y/\ — V v with v — L and R, the basis of the solution is given by two plane-wave solutions 

u(x; e) 
v(x; e) 

and exponentially growing or converging solutions 

u(x; e) 



exp {±ik\\ v x) 




(11) 



v(x; e) 



exp (±k u x) 




(12) 



where 



u„(e) 



y/(l ~ V v f +^ + 1 



v v (e) 



y/(l - V u ) 2 + e^_l_ 
2' 



(13) 



2e T K ' V 2 ^ 

satisfying the normalization |-u^(e)| 2 — |tiy(£)| 2 = 1. The wavenumber k\^ v and the growing 
rate (or converging rate) k v are given by, respectively, 



_L 



^2 [l-V u + y/e 2 + (l-V v y] , 



(14) 



where fcj, 



v/(l-K) 2 + £ 2 -(l-K) 



Wavenumbers fen „ and A;^ which are, respectively, parallel and perpendicular components 
against the normal vector of the wall can be written as k\\ v = k u cos<f) u and kj_ = k u sin<f) u . 
Owing to the translational invariance in y, z directions, we have the law of reflection 



(15) 



Because of the law of reflection, there does not exist retroreflection in the weakly interacting 
condensed Bose system. That is to say, when the Bogoliubov excitation runs toward the 
wall with an incident angle, the reflected excitation does not go back the way one has come. 
We also have a relation between the incident angle and the transmitted angle given by 



k L sin L = k R sin R . 



(16) 



In the low energy limit, the energy has the relation e pa c u k u , where c u is a sound speed of 

Bogoliubov phonon in the dimensionless form given by c u = a/1 — V u . From the relation in 

Eq. (fl6l) . in the low energy limit, we recover the Snell's law 

sin L _ sin 0r 
c l c r 

We assume that an incident wave comes from the left side of the potential barrier as 
shown in Fig. HJ Far from the barrier, exponentially-diverging components e~ KljX and e KRX 
should be absent in the physical solution, and exponentially-converging components e KhX 
and e~ KRX are negligible, so that we have the following asymptotic form: 




\a(e) exp (ik\\ t i,x) + b(e) exp (— ik\\^x)\ 



v h ie 



for x <C — 1 



(18) 



c(e)exp (iku R x) 




for x ^> 1. 



as 



(19) 



We introduce functions S(x,e) and G(x,s) [l, 2, 5, §| 

S(x; e) = u(x; e) + v(x; e) 
G(x; e) = u(x; e) — v(x; e). 

According to Eq. ([TBI) , we have the asymptotic form of the function S(x; e) in the low energy 
limit up to order 0(e 1 ^ 2 ) as 

S(x; e) ~ e- 1 / 2 a/2(1 - V L ) {a(e) + b(e) + ie-f L [o(e) - 6(e)] x} , for 



Sfoe) ~ r 1 / 2 ^!! - V n )c(e) (1 + ie 7R x) 



for x ^> 1, 



(20) 



where we define 7^ = cos 0^/ yl — V^. 

From the Bogoliubov equation in Eq. (flOl) . functions S*(a;; e) and G(x; e) satisfy following 
equations: 

r h + k 2 ± /2\S{x,£) = eG(x,e) 
+ k 2 j2\G(x,e) = eS(x,e), 

where — h + 2\E' 2 (x). We expand functions S(x; e) and G(x; e) with the energy e as 

00 

S{x;e) = e^l 2 £ e n S^(x) 



(21) 



n=0 

00 



G(x;e) = e- 1 ' 2 £ e n G^(x) 

n=0 



(22) 



S 



We note that the factor e l l 2 comes from the normalization in Eq. (fT3"l) . Equations for 
S^(x), G^(x), and S^(x) are given by 

hS^{x) = 
fr-)G( >(x) = (23) 

First, we consider the solution G^. The operator has a second-order ordinary 
differential. The solution of can be written by 



G^(x) = B l9l (x) + B u gn(< 



(24) 



with constants Bi, Bu and two independent solutions g\{x) and gn(x) satisfying, respectively, 

dgi(x) 



01 W 



0, ft (rc = 0) = 1, 



h^gn(x) = 0, (?n(x = 0) = 0, 



dx 
dgu(x] 







z=0 



x=0 



(25) 



At | a; | ^> 1, the asymptotic behavior of can be written as —-^z + 2(1 — V v ). Hence, 
any linear combination of gi(x) and gn(x) generally diverges exponentially at \x\ ^> 1. On 
the other hand, G^ should not have any exponentially diverging terms at \x\ ^> 1. Thus, 
it is allowed that we set B\ = Bu = as in Ref. jg|, and hence we do not have to treat the 
special solution of the inhomogeneous equation of the third equation in Eq. 
We shall consider general solutions of homogeneous equations given by 

hS {1) (x 



hS (0) {x)=0, hS w {x) = 0. (26) 

Equations in Eq. fl26|) are second-order linear differential equations, and hence solutions 
(x) and (x) can be written as 

= A[ 0) Sl (x) + s n (x) 
SV> = Ap ai (x)+A$ aa (x), 

by using two independent solutions si(x) and Sn(x) satisfying 

dsi(x) 



(27) 



hsi(x) = 0, Si(x = 0) = 1, 



hsn(x) = 0, 



sn(x = 0) = 0, 
9 



dx 
dsnix* 



dx 



x=0 



x=0 







1. 



(28) 



At \x\ ^> 1, the asymptotic behavior of h is described as h ~ — \-^, so that we have 



(29) 



solutions as 

s\{x) ~ ai ]V + f3i )V x 
s n (x) ~ an^ + Pn^x, 

with v = L for x <C — 1 and = R for x ^> 1. 

As a result, we have the asymptotic behavior of the solution S(x; e) given by 

+e^ 2 [A 1 (e)(3 ltL , + A u (e)Mx + 0(e 3/2 ), (30) 

with v — L and R, where we define = + eA^ and An(e) = + eA^ and use 

Eq. (J29D. 

Let us use the property that the wave function of the excited state in the low energy limit 
corresponds to the macroscopic wave function of the condensate [7j. Considering solutions 
si,n{x) satisfying hsi^x) = 0, the solution of the Gross-Pitaevskii equation h^/(x) = can 
be written as 

= * «i(aO + %s n (x), (31) 
where = ^(x = 0) and ty' = • Using asymptotic forms for si(x) and sn(x) in 

x=0 

Eq. ( |29i) . we evaluate asymptotic behaviors of the condensate wave function as 

V(x) ~ (* oi^ + *oan, w ) + (tt /V + %Mx, (32) 

with v = \j and R. 

We know that asymptotic behaviors of the condensate wave function are given in Eq. 
([5]), so that we have constraints for coefficients ai„, am-, /3^ u , and /3n v given by 



^oA^ + ^oAi,^ = o, 

where z/ = L and R. The second equation of Eq. ( 1331) shows that vectors t (Pi, u , ftii,v) and 
t (\E'o,^ f o) are orthogonal. Thus, we write the vector t (Pi <u , fln,u) as 

(34) 
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with a constant %. 

Comparing Eq. (|2"0|) and Eq. (1501) . we have equations for coefficients Aj^e), the am- 
plitude transmission coefficient t(e) = c(e)/o(e), and the amplitude reflection coefficient 

r(e) = 6(e)/o(e), 



A I (e)a^ + A n (e)an, L = ^2(1 - V[)a(e)[l + r(e)] + 0(e 2 ) 
A i ( £ )A,l + ^ii(£)Ai > l = ^a(£)[l-r(£)]£CO S L + 0(£ 2 ) 



(35) 



and 



Ai(e)ai, R + ^ n (e)an,R = ^ 2 (! ~ Vk)a(e)t(e) + 0(e 2 ) 
^i(e)A,R + ^ii( £ )Ai,R = iv^a^t^^cos^R + O^ 2 ). 



(36) 



Using Eqs. (133|) . (JMj) . (1351) . and (136|) . we have the amplitude transmission coefficient £(e) 
and the amplitude reflection coefficient r(e) in the low energy limit e — ► 

Hmt(e) = - ^ R - 
V ; Cr + Cl 

(37) 

r / \ Cr - Cl 
^ £ -*° Cr + Cl 

where we define £„ = /?„/ cos^. 

We shall consider the relation between £l and Cr on the basis of constancy of energy flux. 
According to Ref. one has the energy flux averaged over the time in dimensionless form 
given by 

<Q(r)> = |lm [S*(v)VS(v) + G*(r)VG(r)] . (38) 

We shall consider the asymptotic behavior of the energy flux parallel to the normal vector 
of the wall by using Eq. ([TBI . The asymptotic behavior of the energy flux in the x-direction 
for x <C — 1 is given by 

(Q|I, L ) ~ s (\u + , L \ 2 + \v + , L \ 2 ) fc||, L [1 - \ris)\ 2 ] \a(e)\ 2 . (39) 

On the other hand, we have the asymptotic form of the energy flux in the x-direction for 
x ^> 1 given by 

(Qh,r> ~ e (|n +iR | 2 + |S +iR | 2 ) k lltR \t(s)\ 2 \a(s)\ 2 . (40) 
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In the low energy limit, the energy flux can be written as 

(<2||,l) ~ ec L cos(p L [l - \r(e)\ 2 ]\a(e)\ 2 



HI 

|2|„/_M2 



(41) 



<Q||,r> ~ ec R cos(p R \t(e)\ 2 \a(e)\ 2 . 

Regarding the constancy of the energy flux (Q\\ t i) = (Q\\,r) as the constraint for the 
variables Cl and Cr, we obtain the ratio 

Cl = c R cos 0r 
Cr c l cos l ' 

As a result, we have the amplitude transmission coefficient t(e) and the amplitude reflection 
coefficient r(e) in the low energy limit as 

,. s 2c L COS0 L .. C L COS L - Cr COS 0R 

hmt(e) = ■ — , hmr(e) = ; — . (43) 

£^0 C L COS L + Cr COS 0r e^O C L COS (p L + Cr COS </>R 

These amplitude transmission coefficient t and amplitude reflection coefficient r satisfy the 
relation: 1 + r = t, as seen in the classical wave mechanics. The transmission coefficient 
T(s) and the reflection coefficient R(e) are, respectively, given by ratios of the transmitted 
energy flux and of the reflected energy flux to the incident energy flux. In the low energy 
limit, T(e) and R(e) are, respectively, given by 

,. rr>t \ 4c L C R COS L COS 0r 

hmT(e) = — ■ — - 

e^O (C L COS0 L + CrCOS0r) 2 

(44) 

,. n , , (C L COS0 L - CrCOS0r) 2 

\mxR{e) = — ■ — -, 

e^O (c L COS L + Cr COS (p R ) 2 

satisfying the relation: T(e) + R(e) = 1. 

Equation (14411 and the Snell's law in Eq. (1171) describe the transmission and the reflection 
of the Bogoliubov excitation in the low energy limit. We note that T(e — > 0) and R(e — > 0) 
in Eq. (|44l do not depend on the detail of V(x) near x = 0; they depend on V(x) only 
through the asymptotic values Vl and Vr. The perfect transmission T(e — > 0) = 1 and 



R(e -»• 0) = 
earlier studies 



ows 



bl 
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rom cl = cr, or equivalently Vl = Vr. This result is consistent with 
s| . The transmission and reflection coefficients in Eq. (I44p are similar 
to but different from those of the classical sound mechanics, as discussed later. 

The formulation and the result up to the present is adopted to the existence of the 
potential V(x) having the general form. Here, let us compare the transmission coefficient of 
our result with that obtained by the numerical calculation in the presence of a potential V(x). 
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FIG. 3: (Color online) The transmission coefficient T as a function of the energy e, where the 
incident angle (pi, = 0. In the numerical calculation, we assume the potential energy given in Eq. 

The solid and dotted lines in Fig. are the numerical result in the case for Vt, = 3 and 5. The 
dot is our analytical prediction in the low energy limit. 

We shall adopt the potential energy V(x) given in Eq. (jH]). In Fig. [31 we show transmission 
coefficients T in the low energy limit given in Eq. ( 1441) and obtained by numerical calculations 
as a function of the energy e. We set the incident angle as 0l = 0. Solid and dotted lines 
in Fig. [3] are numerical results in the case for = 3 and 5, respectively. The dot is our 
analytical result in the low energy limit. We note that the transmission coefficient obtained 
by numerical calculations reaches our analytical prediction as the incident energy goes to 
zero as seen in Fig. [3j This result means that the transmission and reflection coefficients 
in the low energy limit are independent of the potential barrier. These coefficients depend 
only on asymptotic values of the potential energy. 

III. DISCUSSION 

In this section, we shall discuss the property of transmission and reflection of the Bo- 
goliubov phonon. First, we shall show the phenomena of the Bogoliubov phonon having 
something in common with the electromagnetic wave. These are existences of the Brew- 
ster's angle and the total internal reflection. Second, we shall give an interpretation of the 
anomalous tunneling in terms of the impedance matching. Third, we show the relation 
between the weakly interacting three-dimensional Bose gas and the one- dimensional Bose 
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gas treated as the Tomonaga-Luttinger liquids. We also discuss the Andreev- like reflection, 
critically. Finally, we shall propose future problems on the experimental and theoretical 
sides. 



A. Brewster's Angle and Total Internal Reflection 

In addition to the "ordinary" anomalous tunneling discussed in Refs. |2|, [3|, |8j , we find 
another condition where we do have the perfect transmission of the Bogoliubov excitation 
in the low energy limit. When the incident Bogoliubov phonon has the incident angle </>l,b 
defined by 



tan< 



>L,B 



cl/cr, 



(45) 



we obtain the perfect transmission: T = 1 and R = 0. This incident angle 0l,b corresponds 
to the Brewster's angle for the electromagnetic wave, satisfying the relation 0l b + <Pr = 
tt/2 [lj]. 

The transmission coefficient T in the low energy limit as a function of the incident angle 
0l is shown in Fig. |H We use Vl = and Vr = 0.8 as in Eq. (jH]). The incident angle where 
the perfect transmission occurs is seen in Fig. HI as discussed above. 




FIG. 4: The transmission coefficient in the low energy limit as a function of the incident angle. 
We assume Vl = and Vr = 0.8, consistent with the potential energy given in Eq. ©. 
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Moreover, we also find that the Bogoliubov excitation experiences the total internal re- 
flection when the incident angle satisfies the condition 0l > <Pl,c, where the critical angle 
0l, c is defined by 

sin0 LjC = cl/c r , (46) 

where c L < c R , i.e., V L > Vr. 

In Fig. El the transmission coefficient T in the low energy limit in Eq. (j4"4"]) is plotted in 
the Vr - 0l plane. We assume Vl = 0. The solid line at Vr = represents the line where the 
"ordinary" anomalous tunneling occurs, i.e., the perfect transmission occurs independently 
of the incident angle 0l- On the other hand, the dotted line represents the line where the 
perfect transmission occurs in the same condition as the Brewster's law. As discussed above, 
there is a region where the total internal reflection occurs at Vr < Vl = 0, as seen in Fig. El 




FIG. 5: (Color online) The transmission coefficient in the low energy limit given by our analytical 
result is plotted in the Vr - 0l plane. We assume Vl = 0. The solid line at Vr = represents the 
line where the "ordinary" anomalous tunneling occurs. The dotted line represents the line where 
the perfect transmission occurs in the same condition as the Brewster's law. There is a region 
where the total internal reflection occurs. 

B. Impedance Matching of Bogoliubov Phonon 

The Brewster's law and the total internal reflection of Bogoliubov excitations are naturally 
understood by recalling that expressions in Eq. (jUJ) for the transmission and reflection 
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coefficients have the same forms as those of the electromagnetic wave [15| . The transmission 
and reflection coefficients of an electromagnetic wave going from the medium L to R are 
given by 

_ 4Z L ZrCOS L COS 0r 

(Zr COS L + Z L COS 0r) 2 

(47) 

_ (Z R COS L - Z L COS R ) 2 
(Zr COS L + Z L COS 0r) 2 ' 

when the electric field of the incident electromagnetic wave is perpendicular to the incident 
plane. The symbols Zl and Zr denote the impedance of the electromagnetic wave Z el = 



■sj e'/ /? for each medium. Here e' is the permittivity, and p! is the magnetic permeability. 
When the magnetic permeability for two mediums are equal, Zr/Zl becomes equal to cl/cr, 
where cl and cr are light speeds in left and right mediums. As a result, Eq. (E7J) coincides 
with Eq. (jUJ. We note that the anomalous tunneling discussed in Refs. [l, 2, 3, 8] can be 
then simply regarded as an impedance matching between two identical condensates, where 
the impedance of the Bogoliubov phonon is inversely proportional to the speed of Bogoliubov 
phonon. 

It is worth mentioning that results in Eq. ( )44l) are different from the transmission and 
reflection coefficients of the classical sound wave at the interface between two mediums; at 
the interface, the transmission and reflection coefficients 

4Zg c Z R c cos L cos 4> R 

(Z^ COS (f) L + Z£ c COS 0r) 2 

(48) 

= (Z^COS^l - Z£ c COS0r) 2 

(Z$£ cos (f) L + Z£ c cos 0r) 2 
of the classical sound are given by the same form as Eq. ( 1471) . Here, Z£ c and Z^ are the 
specific acoustic impedance Z ac for each medium. Z ac is defined as the ratio of the sound 
pressure to the particle velocity field and is given by Z ac = pc with a mass density p and 
a sound velocity c. The specific acoustic impedance for Bogoliubov excitations are given 
similarly from the calculation up to first order with respect to u and v. The resultant specific 
acoustic impedance for Bogoliubov excitations is given by the product of the mass density 
of the condensate (oc 1 — V u ) and the sound speed (oc |1 — K,| 1//2 ) for each condensate v =L 
and R. Considering c v oc |1 — V u \ l l 2 in Eq. (|4"4"|) . it is obvious that Eqs. (pf4l) and (j4"8"j) are 
different. 
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The difference of the impedance between the classical sound and the Bogoliubov phonon 
comes from the difference of the boundary condition; in the classical sound wave, which 
is the first sound, say, the hydrodynamic collective mode, the transmission coefficient is 
derived from continuous conditions at the interface for the velocity perpendicular to the 
interface and for the pressure [ly]. On the other hand, in the Bogoliubov phonon, which 
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19] , our result is derived 



is the zeroth sound, say, the collisionless collective mode 
from the constancy of the energy flux, and the property that the wave function of the 
Bogoliubov excitation in the low energy limit coincides with the macroscopic wave function 
of the condensate. 



C. Comparison with Results on Tomonaga-Luttinger Liquids 

We shall discuss the relation between transmission and reflection coefficients obtained 
in our system and those obtained in one-dimensional boson systems. As mentioned in the 
introduction, the perfect transmission in the Tomonaga-Luttinger liquids was investigated 
extensively and intensively in the study of quantum wire. Using the renormalization-group, 
Kane and Fisher showed that the barrier is an irrelevant perturbation, say, the perfect 
transmission occurs, when the Luttinger liquid parameter K satisfies K > 1 [9j. In the 
fermion case, K > 1 corresponds to an attractive fermion system, K = 1 a free fermion 
system, and 1 > K > a repulsive fermion system [^J. Safi and Schulz investigated 
the transport through a one- dimensional wire of interacting electrons connected to leads 



sufficiently longer than the wire [lOj. When one considers the multiple reflections, the 
transmission and reflection coefficients are described by Luttinger liquid parameters in both 
leads, say, these coefficients do not depend on the Luttinger liquid parameter of the wire |lo| . 
Safi and Schulz concluded that the perfect transmission occurs when the Luttinger liquid 
parameters and K R of the left and right leads are the same i.e., = K R (> 1). 

Let us compare the transmission coefficient r and the reflection coefficient 7 for the 
density fluctuation in the system treated in this paper with those in the one-dimensional 
system. These coefficients are defined by the ratio of density fluctuations of transmit- 
ted and reflected waves to the incident wave. In the Tomonaga-Luttinger liquids, the 
transmission and reflection coefficients r and 7 with the multiple reflections, are given by 
r = 2U L K R /[U R (K h + K n )), and 7 = (# L - K R )/(K L + K R ). Here, U L and U R are speeds of 
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collective excitations in the left lead and the right lead. In the repulsive boson case, the Lut- 



tinger 

gas 



iquid parameter K is larger than unity. K — 1 corresponds to the Tonks- Girardeau 



211 ] . In the weakly interacting boson system with a short-range interaction, the Lut- 



tinger liquid parameter is given by K v = ttHaJ p v /mg, where m is a mass, g is a coupling 
constant of the short-range interaction, and p v is the ground state density with v = L for 
the left lead and with v = R for the right lead 22j. The speed of the collective excitation 
U v in the weakly interacting boson system corresponds to that of the Bogoliubov phonon 
c v , say, U u = c v = ^Jgp u /m 23|. When the difference of Luttinger liquid parameters in 
both leads is caused by a potential step, the ground state density p v alone in both leads are 
different. As a result, the transmission and reflection coefficients r and 7 with the multiple 
reflections are given by r = 2cl/(cl + cr) and 7 = (cl — cr)/(cl + cr). 

We shall consider the case for the normal incidence of the Bogoliubov phonon up to first 
order with respect to u and v, on the basis of our results in this paper. The transmission 
and reflection coefficients r and 7 defined as ratios of amplitudes of the density fluctuation 
are given by r = t and 7 = r in the low energy limit, t and r are amplitude transmission 
and reflection coefficients in the low energy limit given in Eq. ( 1431) . As a result, we have 
expressions r = 2cl/(cl + cr), and 7 = (cl — cr)/(cl + cr), at normal incidence 0l = 0. 
We note that this result agrees with the result given by the Tomonaga-Luttinger liquids as 
shown above. 

In connection with the Tomonaga-Luttinger liquids, we shall comment on the difference 
between the reflection of the Bogoliubov excitation and the Andreev reflection. Dynamics 
of one- dimensional Bose liquids has been investi gat ed in several papers, where the negative 
density reflection has been found to occur [24|, [25(. Tokuno et al. studied the dynamics 



of one-dimensional Bose liquids 



241 ] . In the Y-shaped potential, the incident packet splits 



into two transmitted packets running in two branches and one reflected packet running 
back in the branch. They found the negative density reflection. Daley et. al. studied the 
wave packet dynamics, in the one-dimensional optical lattice, propagating across a bound- 
ary in the interaction strength, using the time-dependent density matrix renormalization 



group method 



251 ] . They also discussed the negative density reflection. In these papers, the 



negative density reflection is called the Andreev- like reflection, in the sense that the nega- 
tive density reflection is analogous to the reflection of hole-like excitations at the interface 
between super-normal conductors. 
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This Andreev- like reflection should be however distinguished from the Andreev reflection 
in the following sense. If the condensed Bose system has the Andreev reflection, there 
should exist the retroreflection; when the Bogoliubov excitation runs toward the wall with 
an incident angle, the reflected excitation would go back the way one has come. In the 
weakly interacting Bose system separated by the potential wall, however, there exists no 
retroreflection as shown in Sec. II, and hence there exists no Andreev reflection. The 
reflection is just ordinary and obeys the law of reflection in Eq. (fT5l) . This observation comes 
from the study of the tunneling problem between different densities and at oblique incidence 
in the weakly interacting three-dimensional Bose system. As a result, the Andreev-/i£;e 
reflection in Bose system in the one-dimensional Bose system is different from the Andreev 
reflection. 



D. Future Problems 



The realization of the Bose-Einstein condensation in atomic gases 26|, |27j] has allowed 
us to investigate properties of Bogoliubov excitations experimentally. Using nondestructive 
phase-contrast imaging, it was found that the speed of sound induced by modifying the trap 
potential was consistent with the Bogoliubov theory [28] . Using the Bragg scattering, the 
static structure factor of a condensed Bose gas was measured in the phonon regime [29(. It 
was also reported that the excitation spectrum of a condensed Bose gas agrees with the Bo- 
goliubov spectrum using the Bragg scattering [30]. On the other hand, the transmission and 
the reflection of Bogoliubov excitations are one of the issues that has not been investigated 
yet experimentally. 

It was discussed how to observe the anomalous tunneling in Ref . jsj . If one investigates 
the reflection and the refraction of Bogoliubov excitations, a potential step is needed. The 

n 

potential step could be made by using a detuned laser beam shined over a razor edge |14| . 
Bogoliubov excitation with a sole wave number k could be produced by making use of the 
Bragg pulse {29], Q • An advantage of making use of the Bragg pulse to make Bogoliubov 
excitations is that one can produce an excitation with a sole wave number k having an 
arbitrary incident angle against the potential barrier. 

By making use of the local modification of the trap potential, an excitation can be also 
created. In this case, a density modification is composed of excitations with many modes, 
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and hence the transmission and reflection coefficients should be estimated using a mode 
decomposition, in order to study the energy dependence of the transmission coefficient. 

If we use a box trap [31], these experiments are analogous to those of the reflection and 
the refraction of the wave at the interface between two different homogeneous mediums. 
When we use a harmonic trap, not a box trap, the problem is extended to a kind of the 
problem of the reflection and the refraction with a refractive index which changes spatially. 

In the weakly interacting one-dimensional Bose system with a short-range interaction, it 
is known that the excitation spectrum agrees with that of the Bogoliubov excitation [23[]. In 
this one-dimensional system, we show that transmission and reflection coefficients also agree 
with those of Bogoliubov phonon in the long wave length limit, in the present paper. The 
result on the Tomonaga-Luttinger liquid suggest to us that the perfect transmission occurs 
in the symmetric system even if the interaction is strong. From this suggestion, not only in 
the weakly interacting three-dimensional Bose system but also in the strongly interacting 
three-dimensional Bose system beyond the mean- field treatment, the tunneling problem of 
the collective excitation should be then investigated theoretically and experimentally. The 
strongly interacting Bose gas can be made using the Feshbach resonance 32(]. It is a problem 
whether the anomalous tunneling could be observed or not using experimental techniques 
mentioned above, in such a dilute ultracold gas under the control of the interaction strength. 
On the other hand, it is also a problem whether the anomalous tunneling could be observed 
or not in superfluid He-4 which is a high dense system and a strongly interacting Bose liquid. 

Before closing this section, we shall propose a problem on the theoretical side. Recently, 
Tsuchiya and Ohashi studied tunneling properties of Bogoliubov phonons taking notice of 
the quasi-particle current near the potential barrier 33j . They found that the quasi-particle 
current increases near the potential barrier inducing the supercurrent counterflow. In order 
to conserve the total current, they use the Gross-Pitaevskii equation added in the anomalous 
average. This formulation brings the gapful excitation. Within their formulation with the 
gapful excitation, they did not confirm whether the anomalous tunneling occurs or not. On 
the other hand, it is reasoned that the fact that the wave function of the excited state in the 
low energy limit corresponds to the condensate wave function, say, the gapless excitation 
declared by the Hugenholtz-Pines theorem [34( , is necessary to the occurrence of the anoma- 
lous tunneling. On the basis of this idea, Kato et al. used the Popov approximation in order 
to show the occurrence of the anomalous tunneling even at finite temperatures 8j. However, 
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this treatment is not also necessarily sufficient at finite temperatures, because the Popov 
approximation does not satisfy the conservation law. Hohenberg and Martin [35J showed 



that the Ward-Takahashi relation 
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371 ] which guarantees the conservation law derives the 



Hugenholtz-Pines theorem. Within the theory satisfying the number conservation and the 



gapless excitation, for instance using formulations 38|, |39j, problems whether the anomalous 



tunneling occurs or not and how the quasi-particle current behaves still remain. 



IV. CONCLUSION 



We find that the Bogoliubov phonon experiences the perfect transmission distinguished 
from the "ordinary" anomalous tunneling, when the incident angle satisfies the specific 
condition equal to the Brewster's angle, i.e., the sum of the incident angle and the refracted 
angle is tc/2. We also find that the Bogoliubov excitation experiences the total internal 
reflection. Introducing the impedance for the Bogoliubov phonon, the anomalous tunneling 
can be regarded as the impedance matching. In the weakly interacting Bose system with the 
short-range interaction, the transmission and reflection coefficients are consistent between 
the Bogoliubov theory and the Tomonaga-Luttinger liquid. The negative density reflection 
in the interacting condensed Bose system cannot be necessarily identified with the Andreev 
reflection. 
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